Disentanglement, Bell-nonlocality violation and teleportation capacity of the decaying 

tripartite states 
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Dynamics of disentanglement as measured by the tripartite negativity and Bell nonlocality as 
measured by the extent of violation of the multipartite Bell-type inequalities are investigated in this 
work. It is shown definitively that for the initial three-qubit Greenberger-Horne-Zeilinger (GHZ) or 
W class state preparation the Bell nonlocality suffers sudden death under the influence of thermal 
reservoirs. Moreover, all the Bell-nonlocal states are useful for nonclassical teleportation, while 
there are entangled states that do not violate any Bell-type inequalities, but still yield nonclassical 
teleportation fidelity. 
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I. INTRODUCTION 

Entanglement and Bell nonlocality are two basic ingre- 
dients of a quantum state which are intimately related to 
each other [lj . Their relations have been the research in- 
terests for many years and remain still an important sub- 
ject deserving to be investigated. Entanglement, which 
describes correlations between two or more subsystems, 
refers to the state of a composite system that cannot be 
written as products of states of each subsystem. As a 
physical resource, entanglement is crucial for nearly all 
applications related to quantum information processing 
(QIP). One such application is quantum teleportation 
by which an unknown state can be transmitted from the 
sender to a distant receiver with the help of local opera- 
tions and classical communication (LOCC). However, not 
all the states that are entangled can be used for telepor- 
tation with fidelity (see sections below) better than that 
achievable via classical communication alone, and the fi- 
delity is even not a monotonic function of the degree of 
entanglement of the resource [H-H)]. This demonstrates 
that entanglement may only reveals certain aspects of a 
quantum state. 

Bell nonlocality corresponds to another quantum cor- 
relation of a quantum state that cannot be reproduced by 
any classical local hidden variable models. This nonlo- 
cal property are manifested unambiguously by violation 
of different Bell- type inequalities, which plays a funda- 
mental role in better understanding of the subtle aspects 
of quantum mechanics [l|, @ . Historically, the violation 
of Bell inequalities has been considered as a means of 
determining whether there is entanglement between two 
qubits, for the inseparability of a bipartite pure state 
corresponds to the violation of the Bell inequality in 
the Clauser-Horne-Shimony-Holt (CHSH) form, and vice 
versa @. But this is not the case for the mixed states 
(which are in practice the ones always encountered). As 
demonstrated initially by Werner Q, there exist bipar- 
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tite mixed states which are entangled but do not violate 
any Bell-type inequalities. Further studies also showed 
that the maximal violation of a Bell inequality does not 
behave monotonously under LOCC 0, Q- 

The ability for transmitting information reveals also 
an important aspect of nonseparability of a quantum 
state which is intimately related to entanglement and 
Bell- nonlocal correlations Q. For bipartite two-qubit 
systems, it has been demonstrated that all states that 
violate the CHSH form of Bell inequality can be used for 
teleporting an arbitrary one-qubit state with nonclassi- 
cal fidelity [4J, while there are entangled mixed states 
which do not violate any Bell-type inequalities, but still 
yield nonclassical teleportation fidelity [3j. For the tri- 
partite systems, as we know, there are two distinct classes 
of multipartite entangled states, i.e., the Greenberger- 
Horne-Zeilinger (GHZ) class and the W class, which 
bear incompatible multipartite correlations in the sense 
that they cannot be transformed into each other un- 
der stochastic local operations and classical communi- 
cation @. Particularly, it has been demonstrated that 
besides the Bell states, the three-qubit GHZ and W class 
states in the form of |^ghz) = (|0Q0) + \lll))/y/2 and 
\4) W ) = (\/2|001) + |010) + |100))/2 can also be adopted 
as quantum channels for perfect teleportation under ideal 
circumstance PIE!- Since they both are unavoidable 
disturbed by the surrounding environments, it is natural 
to ask for their robustness against decoherence in terms 
of their teleportation capacity, and its possible relations 
with the degrees of entanglement and Bell-nonlocality vi- 
olation. 

From an applicative point of view, one may hopes that 
quantum correlations which are crucial for QIP can be 
maintained for sufficiently long times to permit designed 
tasks to be fulfilled. But in practice, every quantum 
system is open and susceptible to the unavoidable in- 
teraction with its surroundings [l2r{l4j ]. This may leads 
to decoherence and destruction of correlations. Particu- 
larly, under certain circumstances, the entanglement of 
a bipartite state can even terminate abruptly in a finite 
time, a phenomenon termed entanglement sudden death 
(ESD) by Yu and Eberly [l3| and has been recently con- 
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firmed experimentally |15| . Moreover, due to its prac- 
tical applications in QIP 0, [13], the research interest 
in the decay dynamics of Bell-type correlations has been 
renewed. Analogous to ESD, it has also been demon- 
strated that an initial nonlocal quantum state may loses 
its nonlocal property in a finite timescale under the in- 
fluence of external forces. This phenomenon has been 
named Bell-nonlocality sudden death (BNSD) [H Qj| 
and many theoretical efforts have been devoted to it with 
different contexts (20l - l25j . For instance, Ann and Jaeger 
have demonstrated the occurrence of BNSD for the ini- 
tial three-qubit W class state as measured by its violation 
of the MABK inequality [l8j], as well as for the initial 
generic class of tripartite state as measured by its viola- 
tion of the Svetlichny and WWZB inequalities [l9| . 



II. THEORETICAL FRAMEWORK 

In this paper, we would like to investigate the phe- 
nomena of both ESD and BNSD for tripartite states, 
and their relation with fidelity of quantum tclcporta- 
tion. Our system consists of three qubits (here labeled as 
K = A, B,C), each embedded in a thermal reservoir. To 
focus exclusively on the disappearance of entanglement 
and Bell nonlocality as they arise from the influence of 
thermal noise, we assume the qubits are separated by spa- 
tial distances large enough and thus there are no direct 
interactions between them, that is, every qubit interacts 
only, and independently with its own environment. Then 
under the condition of Markovian approximation, the re- 
duced dynamics of the system state p can be described 
by a general master equation of Lindblad form 12 1 
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where (K = A, B, C) are the damping rates of the 
qubits due to their coupling to the reservoir. The gener- 
ators of decoherence are given by Ck, i = \Jn + loj^ and 
= ^f^cr K , with ° 'k being the raising and lowering 
operators. Cka and Ck,2 describe, respectively, decay 
and excitation processes, with rates which depend on the 
temperature, here parameterized by the average thermal 
photons ft in the reservoir. For the limiting case of van- 
ishing temperature (i.e., ft = 0), only the spontaneous 
decay term survives, leading to a purely dissipative pro- 
cess which drives all initial states to an unique asymptotic 
pure state, in which the three qubits are in their ground 
states. For the opposite case of infinite temperature (i.e., 
ft — > oo), decay and excitation occur at exactly the same 
rate, and the noise induced by the transitions between the 
two levels brings the system into a stationary, maximally 
mixed state. Physically, the above models can describe, 
e.g., three two-level atoms with interatomic separations 
larger than the spatial correlation length of the reservoir 



such that the collective damping and the collective shift 
of the atomic levels are negligible [l2T |. 

Solutions of the above master equation with arbitrary 
initial conditions can be derived exactly in several dif- 
ferent ways I26l. |27H . and here we use the operator-sum 
representation [26[ . In many situations of physical inter- 
est, this representation allows a transparent analysis of 
system dynamics without invoking the explicit forms of 
the initial conditions. In the standard basis expanded 
by the eigenstates of the product Pauli spin operator 
a A ® a B ® a C > tn e reduced density matrix for the three 
qubits together is given by the following completely pos- 
itive and trace preserving (CPTP) map [26| 



P (t) = £{p(0)} = Gijkp(0)G\ : 



ijk> 



(2) 



i,i,fc=l 



where the time-dependent tensor-product superoperator 
C = Ca ® Cb ® C-c contains 64 terms. The Kraus op- 
erator Gijk describes the interaction of the qubits with 
the thermal reservoir and satisfies the CPTP relation 
J2ijk G\jkGijk — 1 for all t. Since the thermal noise 
operates locally on individual subsystems, Gijk can be 
expressed in terms of the tensor products of Ef, Ef 



and Eu as G 
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Here Ef, Ef and 



E^ are the Kraus operators describing time evolution of 
each qubit alone, and individually satisfy the usual com- 
pleteness condition for the operator-sum decomposition 
of CPTP maps 26]. Their explicit forms are as follows 



E* = 
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where the time-dependent factors pk (K — A, B, C) 
appearing in the above equations are given by px = 
e -(2fi+i) 7K t/2. Note that when ft = 0, only Ef and Ef 
survive, so solution (2) has simple analytical form. We 
will apply solution (2) to analyze the effects of the ther- 
mal reservoir on the decay dynamics of entanglement and 
Bell-inequality violation for the initial GHZ and W class 
states. For simplicity, we will take the noise properties 
to be the same for the three qubits such that "/a — 1b — 
7c = 7 and PA = Pb = Pc = P = e~ { - 2n+1 ^ t / 2 . 

After obtaining explicit forms of the reduced density 
matrix p(t), we can discuss dynamics of disentanglement, 
Bell-nonlocality violation and the ability of p(t) for quan- 
tum teleportation. To describe the disentanglement pro- 
cess, we need a concrete measure of entanglement con- 
tained in a quantum state. The tripartite negativity N, 
which was introduced by Sabin and Garcia- Alcaine [28[ 
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is particularly convenient for the case of current interest. 
It can be calculated explicitly from the density matrix 
p(t) as 



classes: 



N = [Na-bcNb-caNc-ab 



,1/3 



(4) 



where N A _ BC = -J^i^f, N B _ C A = -J2i^F and 
Nc-ab = — Yli m£ are the negativities introduced by Vi- 
dal and Werner [29[ , and the sums are taken over all the 
negative eigenvalues pf of the partial transpose p TK (t) of 
p(t) with respect to the subsystem K . Note that for the 
mixed states, N is not able to quantify multipartite en- 
tanglement fully, but its positivity ensures that the state 
under consideration is not separable [28| . Thus if the Bell 
nonlocality of a quantum state dies out before that of the 
tripartite negativity, one can say that it dies out before 
that of entanglement. 

Moreover, we will consider Bell-type nonlocality in dif- 
ferent contexts and we use two classes of multipartite 
Bell-type inequalities to detect the existence of nonlo- 
cal correlations as measured by the extent of their viola- 
tions. The first one we are interested in is the Svetlichny 
inequality |(<S) p (t)| ^5 4 [3(|, which distinguishes gen- 
uinely tripartite Bell nonlocality associated with p(t). 
Here (S) p (t) = tr[«S/?(t)] is the expectation value of the 
Svetlichny operator given by 



S = M A M B M C + M A M B M' C + M A M' B M C 
+M' A M B M C - M' A M' B M' C - M' A M B M C 
-M' A M B M' C - M A M' B M' C , (5) 



where the measurement operators Mk and M' K corre- 
spond to the measurements on each of the subsystems K, 
while the primed and unprimed terms denote the two dif- 
ferent directions in which the corresponding party mea- 
sures (the same applies also to the WWZB operators). 
A quantum state p(t) violates the Svetlichny inequality 
whenever |(<S) p (t)| > 4, and in quantum mechanics the 
Svetlichny inequality is violated up to \(S) p ^\ — 4\/2, 
which is achieved only when the system is prepared in 
the maximally entangled GHZ state @. 

The second inequality convenient for our purpose is 
the WWZB Bell-type inequalities [H|32j]. For the three " 
qubit system, there are 256-element set of such inequal- 
ities, which belong to five distinct classes due to some 
basic symmetries, and the behavior of a single class is 
identical to that of all members of that class, i.e., the vi- 
olation of even a single element of each class is sufficient 
for Bell nonlocality of that class !31]. Thus we only need 
to consider one inequality from each of the five distinct 



Bpi = 2M A M B M C , 

B P2 = ~{-M A M B Mc + M A M B M' C + M A M B M C 
+ M' A M B M C + M a M' b M' c + M' A M B M' C 
+ M A M' B M C + M' A M' B M' C ), 
B P3 = [M A (M B + M' B ) + M' A (M B - M' B )]M C , 
£p4 = M A M B (M C + M' c ) ~ M' A M' B (M C - M' c ), 
B P5 = M A M B M' C + M A M' B M C + M' A M B M C 



M' A M' B M' C , 



(6) 



For nonlocal quantum states, there should be at least 
one of the \{B m ) P (t)\ > 2 (I = 1,2,3,4,5), where 
(Bpi)p(t) = tr[Bpip{t)\. For the behavior of a system 
to be describable by a fully local hidden variable model, 
however, all of the WWZB set of inequalities must be 
satisfied jointly. 

To characterize the quality of the teleported output 
state p ou t under the influence of thermal reservoir, we 
calculate the average fidelity (the fidelity F{6, <p) — 
(Vinlpoutlvin) averaged over all pure input states on the 
Bloch sphere) defined as 



F K 



4tt .,„ 



d9sm0F(0,< 



(7) 



where \ip in ) = cos(0/2)|O) + e^sin(0/2)|l) is the in- 
put state needs to be teleported, with ^ ^ it and 
^ <p $J 2tt being the polar and azimuthal angles, re- 
spectively. When the aforementioned initial GHZ class 
state |V>Ghz) = (|000) + |111))/V2 or the initial W 
class state \ip w ) = (V^|001) + |010) + |100>)/2 is used 
as quantum channel, by following the methodology of 
Ref. [lj], the average teleportation fidelity at an ar- 



bitrary time t can be derived readily as i 7 kv(pGHz) = 
(1 + Pghz 4+55+88 + 2p G 8 HZ )/3 and F av (p w ) = (1 + 
p 22 + 33+44+88+35-46 + 2^%)/^ where the abbrevia- 
tions p «±«±™±- = pg( t ) ± p«(t) ± pT[t) ± , . >f with 

Pn(*) (n = GHZ or W) being the elements of the den- 
sity matrix pn(t). Since the noise imposed by the thermal 
reservoir gives rise both to decoherence and to entangle- 
ment losses, unit fidelity cannot be achieved for this case. 



III. ENTANGLEMENT, BELL-INEQUALITY 
VIOLATION AND TELEPORTATION 
DYNAMICS 

Now we begin our discussion about disentanglement 
dynamics and teleportation capacity for the system pre- 
pared initially in the GHZ class state I^ghz)- The den- 
sity matrix at an arbitrary time can be obtained directly 
from Eqs. (2) and (3), and the nonvanishing elements 
are p 1 ^ 8 (t) (i.e., the diagonal elements) and PqhzW- 
Combination of this with Eqs. (4) and (7), one can de- 
rive the tripartite negativity -/V^ghz) and the average 
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FIG. 1: (Color online) Tripartite negativity iV(pGHz) (solid 
curves) and average fidelity -F av (pGHz) (dashed curves) versus 
jt. For every line style, the curves from top to bottom corre- 
spond to the cases of n — 0, n = 0.1, n = 0.2 and n = 0.3. 



fidelity F av (pGHz) analytically as 
1 



N(pghz) = 2 max {°. V 7 " 2 + P e -P}> 

1 p 3 p 4 (l-p 2 ) 2 

^V(PGHZ) - g + y + J + g (2 - + 1)2 - 



(8) 




FIG. 2: (Color online) Critical tripartite negativities 
^Vc(Pghz) (solid lines) and N c (pw) (dashed curves) versus 
n. 



where the corresponding parameters a and /3 appeared 
in the above equations are given by 




_ (l-p 2 )[2n(n + l)(3p 4 -l) + 2p 4 -p 2 ] 
a ~ 2(2n + l) 3 ' 

_ (l-p 2 )[2n(n + l)(p 2 + l)+p 2 } 
P 2(2n + l) 2 ' ( ' 

with the factor p being defined below Eq. (3). 

Plots of the tripartite negativity -/V(pghz) and average 
fidelity -FkvCpGHz) versus the rescaled time jt are dis- 
played in Fig. 1 with different n. One can observe that 
the evolution of iV(pGHz) shows an exponentially decay- 
ing behavior, and it is very sensitive to the variations of 
the reservoir temperature. From Eq. (8) one can see that 
the entanglement measured by the tripartite negativity 
disappears if a 2 + p 6 ^ /3 2 . When n = 0, this simplifies 
to 1 — p 2 1 . Since for all finite values of 7 the factor 
p approaches zero exponentially only in the infinite time 
limit, the sudden death of N(pqhz) does not happen for 
this special case. This behavior is clearly illustrated by 
the black solid curve shown in Fig. 1. When n 7^ 0, 
however, iV(pcHz) behaves very differently. As can be 
seen from Fig. 1, it ceases to exist in a finite timescale 
te which is shortened gradually by increasing the val- 
ues of n. This implies that the devastating effects of the 
thermal reservoir on entanglement of the system becomes 
severe and severe with increasing temperature. Further- 
more, it should be note that for the initial GHZ state 

I^GHZ), ^OgHz) = N A -BC = N B -CA = N C -AB, thus 



FIG. 3: (Color online) Bell-nonlocality violation of pGHz(i). 
Here the black, red, blue and green curves show respec- 
tively dynamics of S(pgbz), B P2 (pghz), B P3 (p G iiz) and 
Bps(Pghz). For every line color, the curves from right to 
left correspond to the cases of n = 0, n = 0.1, n = 0.2 and 
n = 0.3. 

the sudden death of iV(pcHz) indicates the sudden death 
of N A -bc, N B -ca and N C -ab- 

When considering robustness of the initial GHZ state 
as a quantum channel for teleportation, as can be seen 
from the dashed curves shown in Fig. 1, -F a v(PGHz) de- 
cays to the classical limiting value of 2/3 after a critical 
time tt, which decreases with increasing n, and when 
n = one can obtain 7-Pr = ba[(3 + v / 5)/2]. Moreover, 
one can note that the critical time tt is earlier than 
the death time te of tripartite negativity for any fixed 
n. This indicates that not all the three-qubit entangled 
states generated from the initial GHZ class state are use- 
ful for nonclassical teleportation. In fact, a minimum 
tripartite negativity is always necessary for the achieve- 
ment of nonclassical fidelity for the situations considered 
here. As displayed evidently in Fig. 2, the critical tri- 
partite negativity -/V c (pghz) after which the teleportation 
protocol fails to achieve a nonclassical fidelity is increased 
by increasing the reservoir temperature. 

As mentioned before, Bell-inequality violations may 
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act as an indicator of the usefulness of entanglement. 
Let us discuss now this issue by states evolving in time. 
Our aim is to individuate the Svetlichny inequality vi- 
olation regions characterized by |(S) PGH z(t)l > 4 and 
the WWZB inequality violation regions characterized by 
\(Bpi) PGHZ ( t )\ > 2 (I = 1,2,3,4,5), and compare their 
possible relations with disentanglement and average tele- 
portation fidelity. For this purpose, we calculate the ex- 
pectation values of the Svetlichny operator S and the 
WWZB operators Bpi. For the initial GHZ state, the 
measurement operators for the first qubit (i.e., qubit 
A) are defined as Ma = cr y and M' A = a x 19], while 
the measurement operators for the second and the third 
qubits arc defined with respect to the first one by a ro- 
tation 

( M K \ _( cos0 K -sinO K \ ( M A \ 

\ M' K ) ~ \ smO K cos Ox ) \ M' A ) ' [W> 

where Ok (K = B, C) are the rotation angles. Combina- 
tion of this with Eq. (5) one can obtain 

( 5 )pghz(*) = 4p 3 (sin6» sc - cos6» BC ), (H) 

where Obc = $b + Oc, and the same notation will be 
used throughout this paper. Recall that if |(<S)p GHZ (t)l > 
4, the state PGHz(t) is genuinely tripartite Bell nonlo- 
cal [l]|. Since p = e -( 2n + 1 )7*/2 an( j t^e modulus of 
the trigonometric term takes its maximum y2 when- 
ever Obc = — 7r/4 or 37r/4, the maximum expectation 
value of the Svetlichny operator S evolves according 
to |(S) PGHz(t) | = 4 v / 2e- 3 ( 2fi+1 ^ t / 2 , and approaches the 
classical threshold value 4 in a finite timescale rs = 
ln2/[3(2n + 1)7]. As can be seen from the dynam- 
ics of S'(pghz) = max{|(5) PGHZ (i)| - 4,0} presented in 
Fig. 3, the temperature of the thermal reservoir can 
influence S(pqhz) to a great extent, and the time re- 
gions for genuinely tripartite nonlocality shrink with in- 
creasing temperature. Particularly, in the small yt re- 
gion we find |(<S) PGHZ (t)l decays up to quadratic terms 
in time as |(5) PGHz(t) | = 2^2(a 2 - 2a + 2) + 0(t 3 ), 
where a = 3(2n + l)jt/2. The relative difference be- 
tween the actual and the approximate results, is in prac- 
tice negligible when both are above the classical threshold 

K 5 )pGHz(rs)l = 4- 

Next we extend the analysis of Bell-type nonlocal cor- 
relations in tripartite states addressed by violation of the 
WWZB inequalities. The expectation values for the five 
distinct classes of the Bpi operators for the state pghz (t) 
can also be calculated analytically, which are given by 

( s Pi)pGHz(t) = ( s P4>p G Hz(t) = 2p 3 sin 6»sc, 
(£p 2 ) PGHZ (t) = -p 3 (2 sin BC + cos BC ), 
( S P3>PGHz(t) = 2p 3 (sin6» BC - cosObc), 
(£p5)pghz W = ~4p 3 cos0 B c- (12) 

Recalling that p — e~' 2,l+1 ) 7 */ 2 , one sees immediately 
that the state pghz (t) does not violate the WWZB-type 



inequalities with respect to the operators Bp\ and Bp± in 
the full time region. For the remaining three classes, plots 
of Bpi(pghz) = max{|(Spi) PGHz(t )| - 2,0} (I = 2,3,5) 
versus the rescaled time 7* are displayed in Fig. 3 as 
red, blue and green curves with different n. Clearly, 
they show qualitatively similar dynamical behaviors, the 
most pronounced difference is the time regions during 
which the corresponding inequalities are violated. For 
the inequality of form P2, since the trigonometric term 
is strictly bounded by y/E, the maximum expectation 
value of B P2 evolves as |(S P2 ) PGHz(t) | = V5e~ 3 ( 2 " +1 ^/ 2 , 
and approaches the critical value 2 in a short timescale 
T B P2 — ln(5/4)/[3(2n + 1)7], which decreases gradually 
with increasing temperature (cf. the red curves shown 
in Fig. 3). Second, for the inequality of form P3, the 
trigonometric term is strictly limited by \/2, and the 
maximum expectation value of Bpz is {(Bps) PGHZ ^)\ = 
2v / 2~e -3 ( 2 ™ +1 - >7 */ 2 , who becomes smaller than 2 in pre- 
cisely the same timescale as that of the Svetlichny in- 
equality, i.e., tb P3 — Tg . Finally, for the inequality of 
form P5, the maximum \(B P5 ) PGliz{t) \ = 4 e - 3 ( 2 "+ 1 )T t / 2 
decays to its critical value 2 in a finite time tb P5 = 
ln4/[3(2n+ 1)7]. Note that after time tb P5 all the five 
sets of the WWZB inequalities are satisfied, and the state 
PGHz(i) becomes Bell local, i.e., its correlations can be re- 
produced by a classical local model [l| . This implies that 
under the influence of thermal reservoir, the phenomenon 
of BNSD occurs for definitive times for the initial GHZ 
class state preparation. 

Moreover, by comparing Fig. 3 with Fig. 1, one can 
also note that the lifetime for the tripartite negativity 
is longer than the time region during which the telepor- 
tation protocol outperform those of classical ones or the 
lifetime of Bell nonlocality addressed by violation of the 
full set of WWZB inequalities. This indicates that for the 
initial GHZ state preparation, all the Bell-nonlocal states 
are tripartite entangled and can be used for nonclassical 
teleportation, while there are still some tripartite states 
which are entangled but fail to achieve a nonclassical tele- 
portation fidelity. This conclusion is very similar to those 
of the two-qubit X-type states • 

We now turn our attention to an analysis of disen- 
tanglement, Bell nonlocality and teleportation capacity 
for the initial W state preparation. The density matrix 
Pw(t) can also be obtained directly from Eqs. (2) and 
(3), from which the complete analytical forms of the tri- 
partite negativity and the average fidelity can be derived 
exactly, however, we do not list them here explicitly be- 
cause their expressions are quite involved. Instead, we 
plot numerically in Fig. 4 the time behaviors of both 
N(pw) and F av (pw) for the same values of the parame- 
ters. 

The solid curves in Fig. 4 show the numerical results 
for N(pw) versus jt with different n, from which one 
can note a pronounced difference between the dynamical 
behaviors for the zero and nonzero temperature cases. 
Similar to that for the initial GHZ state preparation, 
the tripartite negativity N(pw) decays exponentially and 
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FIG. 4: (Color online) Tripartite negativity N(pw) (solid 
curves) and average fidelity F av (pw) (dashed curves) versus 
ft. For every line style, the curves from top to bottom corre- 
spond to the cases of n = 0, n = 0.1, n — 0.2 and n = 0.3. 



disappears only in the infinite time limit for the zero 
temperature reservoir, thus there is no ESD happens for 
this special case. For the situations of nonzero temper- 
ature reservoir, however, the tripartite negativity termi- 
nate abruptly in a finite time te, which decreases with 
increasing temperature and is shorter than that for the 
initial GHZ state preparation (cf. Fig. 4 and Fig. 1). 
Thus one can say that the entanglement of the initial 
W state is fragile compared with that of the initial GHZ 
state in the sense that it decays in a faster rate. 

From the dashed curves shown in Fig. 4 one can see 
that in the whole temperature region, the teleportation 
protocol loses its quantum advantage over purely clas- 
sical communication in a finite timescale tt, which is 
significantly shorter than that for the initial GHZ state 
preparation. Thus although they both ensure unit fi- 
delity under ideal circumstance [TO, HH i the channel for 
the initial W class state preparation is fragile compared 
with that of the initial GHZ class state preparation. 
Moreover, from Fig. 4 one can also observe that a non- 
vanishing tripartite negativity N c (pw) is always neces- 
sary for F &v (pw) > 2/3, however, as can be seen from 
the blue dashed curve shown in Fig. 2, N c (pw) decreases 
exponentially with increasing temperature and thus be- 
haves in a markedly different way compared to iV c (pGHz)- 



Now we discuss Bell- nonlocal behaviors of pw(t). The 
measurement operators associated with the fir st q ubit of 
the initial Hastate are Ma = o z and M' A = a x [lg], while 
for the second and the third qubits they are still defined 
as those in Eq. (10). Then the expectation value of the 
Svetlichny operator S can be readily obtained as 



(S) 



Pw(t) 



(2p w 
x (sin 



11+44+66+77+46+47+67-23-25-35 _ j\ 



iBC + COStfsc) 



(13) 




FIG. 5: (Color online) Bell-nonlocality violation of pw{t). 
Here the black, blue and green curves show respectively dy- 
namics of S(pw), Bpz(pw) and Bp^(pw)- For every line 
color, the curves from right to left correspond to the cases 
of n = 0, n — 0.1, n — 0.2 and n = 0.3. 



(4+5\/2)p 4 -(2 + 9\/2/2)p 2 + \/2, which decays with time 
and becomes smaller than 4 in a finite timescale 7x3 = 
ln[(16 + 2072)/(4 + 9a/2 + ^274 + 328^2)] ~ 0.00891. 
For the nonzero temperature reservoirs (i.e., n ^ 0), since 
the density matrix elements of pw (t) are so complicated, 
we define S{pw) = max {\(S) p w (t)\ ~ 4,0} and present 
its dynamical behaviors numerically in Fig. 5 as black 
curves. Clearly, S(pw) decays monotonically with in- 
creasing ft and becomes zero in a finite time rg, after 
which the genuinely tripartite Bell-nonlocal correlation 
disappears, rs decreases with increasing temperature, 
and as can be found obviously from Fig. 5, its magni- 
tude is very small (7x3 < 0.00891), i.e., the genuinely 
tripartite Bell nonlocality for the initial W state prepa- 
ration is fragile and the thermal reservoir destroys it in 
a very short timescale. 

The expectation values for the five WWZB operators 
Spi (I = 1, 2, 3, 4, 5) can also be obtained in terms of the 
density matrix elements of pw(t) as 



For the zero temperature reservoir (i.e., n — 0), the 
maximum value of (S) pw u> simplifies to \(S) Plv (t)\ — 



(Bpi) Pw (t) 
{Bp 2 ) pw (t) 

(BP4) pw(t) 
(#P5> pw(t) 



2) cos 9b cos 9c 



9c) 
^1/2) 



(4/9 n +44+66+77 
+4 ( o^~ 47 sinfls sin 6 C , 

n 25+35-46-47 nna n 

+ Pw - 35+46 - 47 sin(9 B - 

x (cos 9 B c - sm9 B c), 
(2^+ 44+66+77 - 35 + 46 - l)x+ 

+4p w - 47 x-, 
(2p w +ii+66+77 

+4p^" 47 sin<9 B sin (9 C , 

fr, 11+44+66+77+46+47+67-23-25-35 

\ L Pw 

x sin 9 B c- 



V- 



-1) 
(14) 
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where x + = \/2sin(6>B + 7r/4)cos#c, X- = \/2sin(#B — 
n/4:)sm9c, and y± = \/2sm(9c ± 7t/4)cos#b in deriv- 
ing the above equations, we have used the hermiticity 
condition of the density matrix. 

We define B m {p w ) = max{|<£ P i) piv _ w | - 2,0} (I = 
1, 2, 3, 4, 5) to evaluate the extent to which the WWZB- 
type inequalities are violated. Still due to the fact that 
the elements of pw(t) are so complicated, it is difficult to 
express Bpi(pw) compactly with respect to the parame- 
ters -ft and h. Thus we resort to numerical calculations. 
The results show that the inequalities of forms PI, P2 and 
P4 are satisfied in the whole time region (i.e., Bpi(pw) = 
for I = 1,2,4). For the remaining two classes of form 
P3 and P5, examples of decays of Bp 3 (pw) and Bp 5 (pw) 
with different fi are presented graphically in Fig. 5 as 
blue and green curves, respectively. It is obvious that 
they decay monotonously with increasing jt and become 
zero after finite times tb P3 and tb P5 (when n = we 

have7r Sp5 = ln[(20 + 8\/2)/(9 + 2V2+ ^169 + 68^2)] - 
0.10785), both of which decrease with increasing n. The 
time tb P5 demonstrates the time of sudden death of all 
species of Bell- nonlocal correlations for pw (t) under the 
influence of thermal noise. Particularly, by comparing 
the present results with those of PGHz(t), one can note 
that tb P5 here is very small, i.e., the Bell nonlocality for 
Pw(t) is very fragile compared with that of pghz(^). 

Furthermore, one can note that the death time for 
Bell nonlocality is much earlier than that for the tripar- 
tite negativity or the critical time after which the tele- 
portation fidelity F av (pw) < 2/3. This reveals several 
common features. First, there exist time regions during 
which the state pw(t) possesses local correlations even 
in correspondence to high values of tripartite negativity. 
Second, only a small fraction of the tripartite entangled 
states are useful for nonclassical teleportation. But all 
the Bell- nonlocal states yield F^ipw) > 2/3. 

Finally, we would like to emphasize that even when one 
or two of the participating qubits can be well preserved, 
e.g., 7a = or 7^4 = 7s = 0, the tripartite negativity 
and the Bell nonlocality still experience sudden death, 
although in a comparatively longer time. This shows 
that the effects of thermal reservoir on entanglement and 
coherence of a qubit is indeed very different. Since they 
change in the similar manner as that of ja = 1b = 7c = 



7 and no other new results can be drawn, we won't have 
any more discussions about them here. 



IV. SUMMARY 

In summary, we have investigated behaviors of dis- 
entanglement and Bell-nonlocality violation for various 
decaying states and analyzed their relations with ca- 
pacity of these states when being used as quantum 
channels for teleportation. Our system consists of 
three qubits prepared initially in the GHZ or W class 
state and coupled to a thermal reservoir, under the 
influence of which irreversible coherence decay will 
be unavoidable. Depending on the temperature of 
the reservoir, the tripartite negativity can reach value 
equal to zero asymptotically (if h = 0) or at a finite 
time (if n ^ 0). But the sudden death of tripartite 
correlations associated with the Svetlichny inequality 
and nonlocal correlations associated with the WWZB 
inequalities are irreversible in the whole temperature 
regions. Moreover, the tripartite negativity and Bell 
nonlocality for pghz(^) are more robust than that for 
Pw(t) m the sense that they survive in a significantly 
longer times under the influence of thermal reservoir. 
Particularly, pcHz(i) gives a wider time region for 
nonclassical teleportation. Finally, by comparing the 
survival time for Bell nonlocality with the time region 
during which F av (p n ) > 2/3 (II = GHZ or IT), we 
showed that all the Bell-nonlocal states considered in 
this work can be used for quantum teleportation, while 
there also exist a family of entangled mixed states which 
do not violate any multipartite Bell-type inequalities, 
but still yield nonclassical teleportation fidelity. 
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